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The axes of gyroscopes experimentally define local non-rotating frames, i.e. the time-evolution 
of axes of inertial frames. But what physical cause governs the time-evolution of gyroscope axes? 
Starting from an unperturbed FRW cosmology with k = we consider cosmological vorticity per- 
turbations (i.e. vector perturbations) at the linear level, and we ask: Will cosmological rotational 
perturbations exactly drag the axes of a gyroscopes relative to the directions of geodesies to galaxies 
in the asymptotic FRW space? Using Cartan's formalism with local orthonormal bases we cast the 
laws of gravitomagnetism into a form showing the close correspondence with the laws of ordinary 
magnetism. Our results, valid for any equation of state for cosmological matter, are: 1) The drag- 
ging of a gyroscope axis by rotational perturbations of matter beyond the H-dot radius (H = Hubble 
constant) is exponentially suppressed. 2) If the perturbation of matter is a homogeneous rotation 
inside a perturbation radius, then exact dragging of the gyroscope axis by the rotational pertur- 
bation is reached exponentially fast as the perturbation radius gets larger than the H-dot radius. 
3) The time-evolution of a gyroscope axis exactly follows a specific average of the matter inside the 
H-dot radius. In this sense Mach's Principle (that axes of local non-rotating frames precisely follow 
some average of the motion of cosmic matter) is a consequence of cosmology with Einstein Gravity. 

PACS numbers: 04.20.-q, 04.25.-g, 98.80.Jk 



I. INTRODUCTION 



The observational fact. — In tests of general relativity 
in the solar system, two type of things are compared. On 
the one hand measurements of the precession of perihelia 
or of gyroscopes' spin axes in Gravity Probe B [| rela- 
tive to distant stars and quasars. On the other hand the 
solutions of Einstein's equations for the solar system in 
asymptotic Minkowski space, which does not contain dis- 
tant stars or quasars explicitely. In this comparison the 
assumption is made (usually not spelled out) that far 
outside the solar system, but much nearer to us than the 
nearest stars, the local nonrotating frame (of the asymp- 
totic Minkowski space in the solution of Einstein's equa- 
tions) is given by distant stars and quasars. This implicit 
assumption is tested to high accuracy by the comparison 
of the observed perihelion precessions with the predic- 
tions of general relativity . This implicit assumption is a 
basic observational fact. This observational fact has been 
called "MachO" e.g. in 0, although Newton already had 
written that a truly non-rotating system (established ex- 
perimentally by Newton's bucket experiment) is approx- 
imately given by the stars Since stars and quasars 
have proper motions, this formulation cannot be exact. 
But the proper motions of quasars relative to the uniform 
Hubble flow are negligible for present-day tests of Gen- 
eral Relativity. — We conclude that the measurements 
of perihelion shifts and the measurement undertaken by 
Gravity Probe B are tests of two things combined, on the 
one hand tests of Einstein's equations in the solar system, 
on the other hand tests of the principle "MachO" . 
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Classical mechanics, special relativity, and general rel- 
ativity for isolated systems in asymptotic Minkowski 
space give no explanation of the observational fact of 
"MachO" , except saying that this is an accident of initial 
conditions. In these theories one could have different ini- 
tial conditions where all stars could be in rotational mo- 
tion around us relative to our gyroscopes. — Within these 
three theories the local non-rotating frames (one aspect 
of inertial frames) can be experimentally determined by 
axes of gyroscopes, and conversely the time-evolution of 
gyroscope axes is dictated by the laws of inertia, i.e. that 
the gyroscope axes cannot rotate with respect to local in- 
ertial axes. Hence in these three theories things are fully 
consistent but circular, and the question remains: What 
physical cause governs the time-evolution of the axes of 
gyroscopes and inertial frames? 

Mach's Principle. — In the 1880's Mach stated 
clearly and forcefully, as an alternative to Newtonian 
physics, the hypothesis that the axes of local non-rotating 
frames (i.e. axes of gyroscopes) in their time-evolution 
are determined by (are exactly dragged by, precisely fol- 
low ) "some average" of the motion of matter in the uni- 
verse. This is what we take as the formulation of Mach's 
principle. — Many alternative formulations of Mach's 
principle have been proposed by other authors later. We 
shall discuss Einstein's proposal jgj in a future paper. 
Many of the alternatives to Mach's original formulation 
have been enumerated and briefly discussed by Bondi and 
Samuel 2J. Quite a number of these alternatives have al- 
most nothing in common with Mach's ideas. 

Gravitomagnetism. — At the time of Mach there was 
no known mechanism, by which matter in the universe 
could influence the motion of gyroscope axes. With Gen- 
eral Relativity came the needed mechanism, gravitomag- 
netism. Thirring in 1918 analyzed the partial drag- 
ging of the axes of inertial frames inside a rotating in- 
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finitely thin spherical shell with uniform surface mass 
density and total mass M. In the weak field approxima- 
tion, Gn(M/R) 8 heii <C 1, he found that inside the shell 
the axes of local inertial systems at all points rotate rel- 
ative to asymptotic Minkowski space with the same pre- 
cession rate Q — /drag^sheil- For the dragging fraction 
fdrag he obtained / dr ag = ^G N (M/ R) shcn < 1. There is 
only a tiny dragging effect, unless (M/i?) s h c ii approaches 
the value for a black hole. To be relevant for Mach's hy- 
pothesis, exact dragging by the masses in the universe, 
not a small influence, one must go to cosmological mod- 
els. 

In this paper we analyze realistic cosmological models 
(as opposed to toy models) with realistic cosmological 
matter (as opposed to the contrived energy-momentum 
tensors discussed in the literature). We start from an 
unperturbed Friedmann-Robertson- Walker (FRW) cos- 
mology with k = 0, we add the most general cosmologi- 
cal vorticity perturbations (i.e. vector perturbations) at 
the linear level, and we ask: Will cosmological rotational 
perturbations exactly drag the axes of any gyroscope rel- 
ative to the directions of geodesies from the gyroscope 
to galaxies in the asymptotic FRW space? In our anal- 
ysis we cast the laws of gravitomagnetism into a form 
showing clearly the close correspondence with the laws 
of ordinary magnetism. This is achieved by using Car- 
tan's formalism with local orthonormal bases (LONBs) 
and fiducial observers (FIDOs). Our results, stated in 
the abstract and presented in sections IV and V, show 
that Mach's Principle (that axes of local non-rotating 
frames precisely follow some average of the motion of cos- 
mic matter) is a consequence of cosmology with Einstein 
Gravity. The crucial equation is Einstein's Ggj-equation, 
Ampere's law for gravitomagnetism, Eq. 1|28|) in a spa- 
tially flat FRW universe. The mathematical statement, 
what average of the energy flow out there in the universe 
determines the time-evolution of gyroscopes' axes here, is 
given in Eq. i|32|l . A short version of this work appeared 
in Q . In a subsequent paper we shall present the analysis 
for vorticity perturbations on a FRW background with 
k = ±1. — We use the conventions of Misner, Thornc, 
and Wheeler [jj. 



II. VORTICITY PERTURBATIONS, 
FIDUCIAL OBSERVERS, 
AND THE GRAVITOMAGNETIC FIELD 

Cosmological Vorticity Perturbations. — Linear cosmo- 
logical perturbations, see J. Bardeen decouple in 
three sectors, 3-scalars (density perturbations), 3- vectors 
(vorticity perturbations), and 3-tensors (gravitational 
waves). In the vector sector all quantities must be con- 
structed from a 3- vector field with vanishing divergence. 
Hence the 3-scalar Sgoo is zero, the lapse function is un- 
perturbed, and the slicing of space-time in slices E t is 
unique, i.e. there is no gauge ambiguity about the time 
coordinate. For linear vector perturbations the intrinsic 



geometry of each E t remains unperturbed. This holds be- 
cause the perturbations of the Ricci scalar, of the space- 
space components and the trace of T M!y , and of the Ricci 
tensor all must vanish, and for d — 3 the Riemann tensor 
can be built from the Ricci tensor. Our choices of spatial 
coordinates are Cartesian for k — resp spherical FRW 
coordinates for k = 0, ±1. The line element is 

ds 2 = ~dt 2 + {ahi) 2 {dx l f + 2(ah t ) 2 /3 i dx i dt, (1) 

where (3 l is the shift 3-vector. Geodesies on E t are 
straight lines on our choice of chart. We consider vec- 
tor perturbations in an asymptotic FRW universe, i.e. 
f3 l — > for r — > oo. Our coordinates are fixed to "distant 
galaxies" , i.e. to galaxies in the asymptotic unperturbed 
FRW space, and the basis vectors in the coordinate ba- 
sis, e;(P) = d/dx l , point along geodesies in E t from P 
to fixed "distant galaxies" . 

Fiducial Observers. — Our aim is to obtain the laws 
of linearized gravito magnetism in a form analogous to 
electromagnetism in a 3 -I- 1 formulation. What is the 
operational definition for E g (gravitoelectric field) and 
B g (gravitomagnetic field)? According to the equivalence 
principle for a free-falling, non-rotating observer there are 
no gravitational forces at his position, E g = 0, B g = 0. 
It all depends on the choice of fiducial observers, FIDOs, 
with their local ortho-normal bases, LONBs, see Thorne 
et a lfll| . Hence we work in the formalism of E. Car- 
tan ||. Our choice of FIDOs: The world lines of our 
FIDOs are at fixed x l in our coordinates, which are fixed 
to distant galaxies in the asymptotic FRW universe, and 
Cq(P) = mfido(-P)- Hats refer to LONBs, and bars desig- 
nate space-time vectors. We choose the spatial basis vec- 
tors of our FIDOs, ej(P), fixed to directions of geodesies 
on Tit from P to distant galaxies in the asymptotic FRW 
universe. Specifically we fix ej(P) in the same spatial di- 
rections as e,(P) = d/dx l , i.e. in 4-space the directions 
of ej and ej differ by a pure Lorentz boost, see Eq. JSJ be- 
low. The 3-velocity of our FIDOs relative to the normals 
on E t is equal to the shift 3-vector /3* in Eq. QJ. 

The operational definitions of E g and B g . — These def- 
initions are independent of perturbation theory. They 
involve FIDOs (of any given choice) measuring the first 
time-derivatives along their world lines, on the one hand 
of the momentum components p\ of free-falling qua- 
sistatic test particles, and on the other hand of the spin 
components S\ of gyroscopes carried along by the FIDOs, 

^jJP'i = rnE? free-falling quasistatic test particle, (2) 

^-S^-^lBgAS}; gyro comoving with FIDO, (3) 

where t is the local time measured by the FIDO. Arrows 
denote 3-vectors in the tangent spaces spanned by the 
spatial legs of our LONBs. E g = g is the gravitational 
acceleration of free-falling quasistatic test particles rela- 
tive to the FIDO. Eqs. (|2I3|) are the same as for a clas- 
sical charged spinning test particle in an electromagnetic 
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field except that q is replaced by m, and the gyromag- 
netic ratio q/(2m) is replaced by 1/2. Eq. © gives the 
angular velocity of precession of the gyroscope's spin axis 
relative to the axes of the FIDO, 
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gyro 



alphabet both for spatial LONBs (with hat) and for spa- 
tial coordinate bases (without hat). The dual bases (ba- 
sis 1-forms) ¥ for LONBs are defined by (# a ,e s ) = Sf, 
where tildes designate space-time 1-forms. The LONB 
1-forms 9 a are expanded in the coordinate basis 1-forms, 
(4) § a = dx a , i.e. 9 & = (9 & ) a 6 a , by 



which is an equivalent operational definition of B g . 



III. CONNECTION COEFFICIENTS AND THE 
EQUATIONS OF MOTION FOR MATTER 
IN GRAVITOMAGNETISM 



Connection 1-forms. — The connection 1-forms (I>°- 

J b 

resp their components in LONBs (= Ricci rotation co- 
efficients), (u a ~)c, are defined by 



(5) 



In words: the Ricci rotation coefficients (w c -)„ give the 
rotation resp the Lorentz boost (^ c g) of the LONBs rel- 
ative to parallel transport along e"a- Parallel transport is 
given by free fall for mfido = eg and gyroscope axes for 
eb. Relative to FIDOs the equations of motion for free- 
falling test particles (geodesic equation) and for spin axes 
of gyroscopes (Fermi transport) specialized to gyroscopes 
carried along by a FIDO are 



dp" 



tdx 



dS l 



-^ + K)a^3T=0, — + K.) ^=0. (6) 



dt 



dt 



dt 



With Eqs. © the operational definitions Eqs. 10 get 
translated into the equivalent definitions involving con- 
nection coefficients with a displacement index 0, namely 
a Lorentz boost lu-^ per unit time (acceleration g = E g ) 
resp a rotation angle wj-- per unit time (angular velocity 
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gyro 



- ~l B g)i 



(7) 



where By = e rfk B k . 



The computation of connection coefficients in Cartan's 
formalism. — The first step is to express our choice of 
LONBs e a (P) in terms of the coordinate bases e a (P) = 
d/dx a , i.e. ea = {e a ) a e a . To first order in ((3 l /c) 



= 6° - (3 k e\ 



~9 k = ah k 9 k . (9) 



The coefficients of the inverse expansion, i.e. coordi- 
nate bases in terms of LONBs, are (e a ) a = (0 a ) a resp. 

(0 a )a = (ea) a - 

The second step is computing the exterior derivative d 
of the basis- 1-forms, (d9 c ) a p = d a {9 c )p — dp(9 c ) a , where 
[a/3] must be in the coordinate basis. Then one converts 
to components [ab] in the LONB, 

(d9 e ) &i ee -C,J = (e & ) a [d a (9 £ ) - d (9%](e h f. (10) 
The coefficients C- c are identical to the commutation 

ab 

coefficients of the basis vectors, [ea,eg] = C^ c e~c- This is 
easily shown by using d a (9 c ,ef) — = [d a (9 c ) p]^) 13 + 

The third step is obtaining the connection coefficients 
from the commutation coefficients. The definition of the 
connection via basis 1-forms is (y a 9 c )p — —(u} c ^) a (9 d )p. 
We take both the displacement index a and the equa- 
tion's component index (3 in the coordinate basis, and 
we antisymmetrize in [a(3]. This makes the Christoffel 
symbols T p a/3 on the left-hand side disappear, since they 
are symmetric in a, [3. Hence the left-hand side reduces 
to (d9 c ) a p. Dropping the equation's component indices 
[a, 13] gives 



d9 c 



-W C ;Al 



(11) 



This is Cartan's first equation. The wedge product (exte- 
rior product) of two 1-forms a and p is (aAp) a /3 = a a p/3 — 
o~/3Pa- Taking Cartan's first equation in LONB compo- 
nents gives the commutation coefficients, Eq. (I1UII . on 
the left-hand side, and the right-hand side simplifies in 



LONB because 



5f. Hence Cartan's first equation 



in LONB components is 

C ai = K)*- - (^)b- (12) 

This equation is easily solved for the rotation coefficients, 



-r-(e* +(3 k e )- (8) 
ah k 



( W cS)« — ~^Pcba + ^cab ~ ^bat\- 



(13) 



For FRW with k = in Cartesian spatial coordinates 
all hi = 1. For FRW in spherical coordinates and 
k = 0,±1 we have h x = 1, h$ = R(x),h,f, = i?(x)sin# 
with R(x) = {Xi sm Xi smn x}- Since the spatial LONBs 
point in the same spatial directions as the spatial coor- 
dinate bases, we use latin letters from the middle of the 



Connection coefficients for vorticity perturbations on a 
Minkowski background. — To first order in [3 l and with 
Cartesian spatial coordinates on E t , the commutation 
coefficients are very simple to compute, because only 
(9°)i — —f3i is space-time dependent, and the prefac- 
tors in Eq. ljTU|) . (e a ) a and (eg)' 3 , can be set to 1. Hence 
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C ° 

ab 



(ef/3) ab , and 



No) 6 
Nj)o 



- E? 
1 



2 l 3 k k 



No 



0. 



(14) 



All components of connection coefficients with respect to 
LONB's are directly measurable (in contrast to Christof- 
fel symbols, which refer to coordinate bases). 

Connection coefficients for vorticity perturbations on 
FRW with k — 0, ±1. — It is again straightforward to 
compute the commutation coefficients and the connection 
coefficients using Eqs. (0 El EH EH) , 



No) 6 

Nj)o 
No)j 
Nj)fc 



-Eg = -j-ftft, 
1 ahi 

2 ^ k k 2a 2 hihj 



(d0)i 



(15) 
(16) 
(17) 
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^(Hp j + d 3 Li) - ^(HPi + diLj), (18) 



ah 



where Li = log hi. Since we work to first order in the 
vorticity perturbations (i.e. in (3 1 ), we can identify E g 
and B g with vectors in S t . From Eqs. Q15I I16fl . we see 
that the shift vector j3 must be identified with the grav- 
itomagnetic vector potential A g . From Eqs. (|T5l HB|) 
follow 
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B g = cmlAg, Eg = --dt(aAg), 
cmlEg + -Ld t {a 2 B g ) = 0. 



(19) 
(20) 



These equations are identical with the hogeneous equa- 
tions for electromagetism in FRW space-times with k = 
0,±1. 

Equation of motion for free-falling test particles. — The 
equation of motion (geodesic equation) for test particles 
of arbitrary velocities v < c in linear vorticity perturba- 
tions on a Minkowski background reads 



e[E g + (vABg)] h 



(21) 



identical with the one for electromagnetism, except that 
the charge q is replaced by the energy e of the test par- 
ticle. With Eq. and in a stationary gravitomagnetic 



field (E g = 0) Eq. HE) becomes ^(pj) 
the Coriolis force law. Note that Cl syio is minus the rota- 
tion velocity of the FIDO relative to the gyroscopes' axes. 
A homogeneous gravitomagnetic field can be transformed 
away completely by going to a rigidly rotating coordi- 
nate system, i.e. physics in a homogeneous gravitomag- 
netic field is equivalent to physics on a merry-go-round in 
Minkowski space. — Note that there are no terms bilinear 
in v for test particles of arbitrary velocities v < c. 



For FRW with k = we obtain 
\jt {aPl) = ^ 9 + (v^B g ) + HvA0Av)} l . (22) 

For FRW with k ± 1 there are additional terms from 
diLj in Eq. i|18|) . These terms are present even in the 
absence of vorticity perturbations and of Hubble expan- 
sion, because in the spherical basis spatial LONBs are 
not parallelized. 



IV. EINSTEIN'S EQUATIONS FOR 
GRAVITOMAGNETISM: AMPERE'S LAW 

Curvature. — The curvature 2-form lZ a ; has LONB 

b 

components (7?. a -),s d -, which are the LONB components of 
the Riemann tensor R a :,~,. The Riemann tensor can be 

bed 

operationally defined by the action of (V 7 V,5 — V,5V 7 ) 
on the LONB 1-form 9 & , 

(v 7 v 5 - v 5 v 7 )<r = -9\n% s , (23) 

where the covariant derivatives V 7 and V5 must be in the 
coordinate basis. To compute the curvature 2-form we 
first use Vs0 a = — (uj a & ) 9 e . Then we let this right-hand 
side be acted on by V 7 . This gives two terms. One term 
comes from V 7 acting on 9 e , and after antisymmetriza- 
tion in [76] it produces — (uj a g Au e ^) 1 $9 b . The other term 
comes from V 7 acting on the expansion coefficient (num- 
ber field) (u> a g)s, where it can be replaced by <9 7 , and af- 
ter antisymmetrization in [715] it produces — (duj a -^) 1 $9 h . 
Hence we obtain 



b b 



(24) 



which is Cartan's second equation. 

Cartan's 2nd equation in LONB. — To obtain the 
LONB components of the first term of the right-hand 
side, (duj a ^)-j, we must first convert the connection com- 
ponents of Eqs. (fT5l - IT%|l from the LONB to the coor- 
dinate basis, (w%)s — O*A)<|(0 <i )*! then take the exterior 

derivative, <9 7 {(w a g) ^{9 d )s} — [7 *-* S], and then convert 
the result back from coordinate components to LONB 
components. The partial derivative of the product gives 
two terms, one with a partial derivative of (uA)j, the 

other with d~ f (9 d )s, which produces another connection 
1-form component. In the second term of Cartan's sec- 
ond equation these conversions from LONB to coordinate 
basis and back again cancel, since there is no derivative in 
between. The result is Cartan's 2nd equation in LONB 
components, 



[c <-> d] . 



(25) 
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Einstein equations for vorticity perturbations of 
Minkowski space. — For linear vorticity perturbations of 
Minkowski space (with Cartesian coordinates for 3-space) 
all non-zero connection coefficients in Eqs. Q15I- 118ft are 
of first order in the perturbations. Therefore the second 
term of Cartan's second equation can be neglected, and 
in the first term one need not distinguish components 
in LONB from components in the coordinate basis. For 
vorticity perturbations the important Einstein equation 
is the equation for G^ — 

= = 0Msj)t3 = \{™t\b) % . (26) 

Hence Einstein's Ggj-equation for vorticity perturbations 
in Minkowski space reads 

curLB 5 = -167tGatJ £ . (27) 

ji = rpo-i _ ^pj^p^yi j g energy current density, which 
is equal to the momentum density. Eq. Q27JI is identi- 
cal to the original law of Ampere for magnetism, except 
that the charge current J q is replaced by the energy cur- 
rent J e , and the prefactor 4-7T is replaced by the prefactor 
(—16itGn)- In contrast to the Ampere-Maxwell equation, 
the Maxwell term [ptE] is absent in gravitomagnetody- 
namics. The Gjg is an equation at fixed time, a constraint 
equation, called momentum constraint, since the momen- 
tum density appears on the right-hand side of Eq. (|27|1 . 

To see the analogous structures of gravitomagnetism 
and electromagnetism, it is more instructive to formulate 
this constraint equation, as we have done in Eq. I|27|) . via 
the connection 1-forms, which involves (di(3j — djPi), i.e. 
the gravitomagnetic field, than via the extrinsic curva- 
ture tensor Kij, which involves {difij + dj[3i). Of course 
the resulting constraint, if written in terms of A g = [3, is 

the same, AA g — WttGnJe- 

The Gqq equation with the source Tqq is trivially ful- 
filled, since these objects are 3-scalars and therefore van- 
ish in the vector sector. — The source T« vanishes, since 
it is of second order in the perturbation. The G-^ equa- 
tions give do(di(3j + djfli) — 0, i.e. the shear of the field 

(3 has vanishing time-derivative. 

Einstein equations for vorticity perturbations of spa- 
tially flat FRW space. — With Cartesian comoving coor- 
dinates for flat 3-space there are two new terms in the 
connection coefficients, (oj^)? rw = S-^-H and (wm)^ = 

HWj-SjfA). Computing % = (ft Qj ) ?j with Eq. E3 
we obtain the corresponding Einstein equation, 

curLB 9 - 4HA g = -1&tiG n J £ , (28) 

where we have used (3 = A g . The scale factor a of the spa- 
tially flat FRW universe does not appear in these equa- 
tions. From H — ~AttGn{p + p) we see that H < 
for p > —p. Therefore we define the H— dot radius by 
R\ = {-H)-\ and we define p 2 = -4H = (l^) -2 . 

We insert the vector potential A g = [3, we use div A g = 



0, hence curl curl A g = —A A g . Therefore Eq. (|2*5)l be- 
comes 

(-A + fi 2 )A g = -16 7tGjv J e ■ (29) 

The new term on the left-hand side, {—AH (3) = (n 2 A g ), 
dominates for superhorizon perturbations. 

V. MACH'S PRINCIPLE 

Our first result. — The solution of Eq. I|29|) is the 
Yukawa potential for A g — (3 in terms of the sources J e 
at the same fixed time, 

l g (r,t) = -4G N [ d 3 r'J e (r',t) ^ti^ll . (30 ) 
J \i — r'\ 

This is analogous to the formula for ordinary magne- 
tostatics except for the exponential cutoff. The Green 
function which is exponentially growing for r' — > oo is 
rejected on the standard grounds of field theory. The 
Yukawa potential in Eq. I|3U|I has an exponential cutoff 
for \ f — r'\ > 1/ fi. This gives our first important conclu- 
sion: The contributions of vorticity perturbations beyond 
the H— dot radius are exponentially suppressed. 

Our second result concerns the exact dragging of gy- 
roscope axes by a homogeneous rotation of cosmological 
matter out to significantly beyond the if— dot radius (for 
the exponential cutoff to be effective) . This holds for any 
equation of state. This is easily seen from Einstein's Gqj 
equation Ij29(l in fc-space for superhorizon perturbations, 
fcp hys -c {—H), where the A-term can be dropped. Us- 
ing J e = (p + p)vfiuid and H = —inG N {p + p) 
we see that all the prefactors cancel, and we obtain 

(3(x) = -u flui d(f). With Ogyroscope = |(v x (3) and 

with f2fl u id = |(V x i/fluid) we obtain the result that for 
-Rpert ^> Rfj there is an exponentially fast approach to 

^gyroscope — ^matter- (^-Q 

This proves exact dragging of gyroscope axes here by 
a homogeneous rotation of cosmological matter out to 
significantly beyond the H radius. 

Our third result concerns the most general vorticity 
perturbation in linear approximation, and it states what 
specific average of energy flow in the universe determines 
the motion of gyroscope axes here at r = 0. We take 
B g = curing of Eq. (|3T7|) . we set r — 0, which gives 
B g (r = 0) = — if2 gyro , and obtain the equation for f2 gyro 
in terms of the sources at the same fixed time, 

^gyro = 

2G N (p + p)J d 3 r 1 [(1 + M r) e~» r ] [r A v(r)} . (32) 

This is the precise expression for Mach 's principle, it says 
exactly what average of the motion of energy in the uni- 
verse determines fL V ro- Mach had asked: "What share 
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has every mass in the determination of direction in the 
law of inertia? No definite answer can be given by our 
experiences" Q. — In the integrand we have the grav- 
itomagnetic moment density ft g = |(F A J e ), analogous 
to the magnetic moment density of an electric current 
distribution. The expression (FA J e ) — (L e ) is the mea- 
sured angular momentum. This is the lowest term, the 
£ = 1 term, in the multipole expansion of the source for 
Jobs = fgyro = and r SO urce > 0. Higher multipoles can- 
not contribute to the gravitomagnetic field B g at r = 0. 
At each radius r only a term equivalent to a rigid rota- 
tion with angular velocity (i(r) contributes in Eq. (|32|l . 
Using this we obtain 

-* 4 f dr r -* 

^gyro = 77 / n 9 ^matter (f) C ^ (r) , (33) 
3 J H 

where C M (r) is the cutoff function, i.e. the first square 

bracket in Eq. (|3"2|) . For the special case flatter (r) inde- 
pendent of r we recover Eq. (I31|l . 

The r-dependence of the weight function of the energy 
current J £ (r) in the integral for B g (0) for rCi!^ is the 
same (1/r 2 ) law as in Ampere's law. Hence the weight 
function for the measured angular momentum (F A J £ ) is 
(1/r 3 ). For fixed CI the energy current density increases 
linearly with r, and the weight function per unit r is 
the first power of r, Eq. 133|) . Therefore a perturbation 
which is a rigid rotation out to i? pcr t and zero outside has 
a dragging fraction /drag, which grows quadratically with 
i? pcrt , until /drag reaches a value near 1 for i? pe rt = Rjj- 
For i?pert increasing beyond i?^ the dragging fraction 
approaches the value 1 exponentially fast. 

Dark energy with pj p — —1, i.e. a cosmological con- 
stant, does not contribute in Mach's principle, Eq. (|32[1 . 
since there is no flow of energy associated with it, its 
energy current J £ — (p + p)v vanishes. 

Measuring everything relative to axes of gyroscopes at 
a given location (instead of relative to galaxies in the 
asymptotic FRW space) makes the the left-hand side of 
Eq. (|32[1 vanish. Hence Eq. (|32[1 reduces to the statement 
that the angular momentum of matter, measured relative 
to the gyroscopes at the given location, will vanish after 
averaging with the weight r~ 3 and with the exponential 
cutoff Cfj,(r). From the point of view of measurements 
it is preferrable to mesure relative to gyros at one given 
location. But to see the structure of gravitomagnetism 
most clearly, it is best to measure relative to galaxies in 
the asymptotic FRW space. 

Analogous dragging effects in magnetostatics (e.g. a 
rotating charged spherical shell acting on magnetic dipole 
moments inside), have the opposite sign in Ampere's law 
compared to gravitomagnetism, Eq. (|27l) , and this causes 
antidragging in the magnetostatic case. 

Einstein's objection to Mach's principle in his autobi- 
ographical notes of 1949. — Einstein wrote [lj: "Mach 
conjectures that inertia would have to depend upon the 
interaction of masses, precisely as was true for Newton's 



other forces, a conception which for a long time I consid- 
ered as in principle the correct one. It presupposes im- 
plicitly, however, that the basic theory should be of the 
general type of Newton's mechanics : masses and their 
interactions as the original concepts. The attempt at 
such a solution does not fit into a consistent field theory, 
as will be immediately recognized." We have shown, how 
this apparent difficulty is resolved in General Relativ- 
ity, specifically in weak Gravitomagnetism: The relevant 
Einstein equation has the form of Ampere's law with a 
Yukawa term, Eq. 128|) , hence the measured mass-energy 
flow out there in the universe does indeed determine the 
precession of gyroscope axes here, Eq. (|3~2)l . 

VI. MEASURED MATTER INPUT FOR 
EINSTEIN'S EQUATIONS AND FOR 
MACH'S PRINCIPLE 

The input on the right-hand side of Eq. which 
expresses Mach's principle, is the measured angular ve- 
locity of matter (stars, galaxies, etc), measured relative 
to galaxies in the asymptotic FRW space ("asymptotic 
galaxies"). No knowledge of the metric perturbation [3 
is needed when determining the input for Eq. I|32|l . This 
is a purely kinetic input, which means that it is directly 
determined by the measured state of motion of matter 
at a given time without knowing the metric perturbation 
(3. Similarly the measured angular momentum is a purely 
kinetic input, (p + p)(f A v)r 2 dr = (p + p)r 2 flr 2 dr. The 
dynamical output of solving Einstein's Ggj equation is 
(3 = A g , the gravitomagnetic vector potential, in Eq. I|29|) 
and the gravitomagnetic force B g in Eq. I|32|l . The mea- 
sured, kinetic angular momentum must be distinguished 
from the canonical angular momentum, which we intro- 
duce (review) in the following paragraphs. 

Action and Lagrangian for Gravitomagnetism. — The 
Einstein-Hilbert action for linear vorticity perturbations 
on a Minkowski background and for point particles is 

S = (16-itGn)- 1 J d 4 x^ (curli*,,) 2 

+ / dt^2[-mx n + mx n A g (x n ,t)] (34) 

n 

Except for the prefactor (167rGAr) -1 and the sign of the 
first term, the action is the same as for electromagnetism 
without the (9(A) 2 -term. — Generally the equations of 
motion are given by the Lagrangian via the standard 
Euler-Lagrange equations (as in classical mechanics and 
classical electrodynamics), if and only if the Lagrangian 
is defined by 

S = J dtL (35) 

without any metric factors in the integrand. Hence the 
Lagrangian for a point particle in a gravitomagnetic field 
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is given by the square bracket of the matter term in 
Eq. J3U). 

The canonical momentum is defined by (p C an)fc = 
dL/dx k , where k = 1,2,3. From Eq. QgJ we obtain in 
Cartesian coordinates 

Pcan = m(x + A g ). (36) 

This is the same equation as in classical mechanics for 
point particles in an electromagnetic field, except that 
the electric charge q is replaced by m. — The kinetic mo- 
mentum is rax. It is directly measured, it can be used 
as the input for solving Einstein's equations, because it 
is independent of the gravitomagnetic vector potential 
A g = fj, which is an output of solving Einstein's equa- 
tions. On the other hand the canonical momentum (for 
a given measured state of motion) depends on the grav- 
itomagnetic vector potential A g . Therefore the canonical 
momentum cannot be used as a input for solving Ein- 
stein's equations. The canonical momentum cannot be 
determined by a FIDO from measurements before having 
solved Einstein's equations. — In curvilinear coordinates 
of 3-space the general definition of the canonical momen- 
tum, (p ca n)fc = dL/dx k , gives 

(Pca„)fc =m gkn {x n + A n g ) k,n = 1,2,3. (37) 

From its definition via the Lagrangian, the canonical 
momentum is a 1-form in 3-space, i.e. it has a lower 
3-index. — On the other hand we can also start from 
the 4- velocity u v , which is the archetype of a 4- vector, 
multiply with the mass to obtain the 4-momentum p v , 
and pull down the 4-index with ^ g Kl/ . This gives pu — 
gkv{mu 1 ') = mgkn(x n + Ag), which is the same as the 
canonical momentum in Eq. 137f) . — Going to spherical 
coordinates we note that Prf, has the physical meaning of 
canonical angular momentum. 

In the special case of azimuthal symmetry, i.e. when 
= £ is a rotational Killing vector, the canonical angu- 
lar momentum p^ = (p, £) is conserved, while the mea- 
sured, kinetic angular momentum of matter, (rp^), is not 
conserved, because of the gravitoelectric induction field 
E g , Eq. I|2(J|) . which acts in the (^-direction. The canoni- 
cal angular momentum is relevant for the time-evolution, 
i.e. for the dynamics, and conservation laws, not for kine- 
matics (measurements at a given time). 

In the continuum description of matter (energy cur- 
rent density) the LONB components T^ can be used as 
input for solving Einstein's equations, since they can be 
measured without knowing the output /3r of Einstein's 
equations. On the other hand the coordinate-basis com- 
ponents T° k cannot be used as an input, because they 
cannot be determined by measurements without a knowl- 
edge of [3 k - 

In Mach 's principle the input is the observations of the 
angular velocities of stars and galaxies, and from there 
the measured kinetic angular momentum (p + p) [f A v] 
as shown in Eq. I|32|l . — Bicak et al 01 have proposed 
to use p<j,, the canonical angular momentum as the input 



on the right-hand side of the Einstein equation. For the 
reasons given above in this section, we consider this to be 
a fundamental mistake. If one does this, the H A g — term 
on the left-hand side of the Einstein equation (|28() gets 
cancelled, and the exponential suppression factor disap- 
pears. 

Einstein's objection to Mach's Principle in his letter 
to Felix Pirani of 2 February 1954 which is quoted 
by Ehlers in [l^: "If you have a tensor T^ and not a 
metric, then this does not meaningfully describe matter. 
There is no theory of physics so far, which can describe 
matter without already the metric as an ingredient of 
the description of matter. Therefore within existing the- 
ories the statement that the matter by itself determines 
the metric is neither wrong nor false, but it is mean- 
ingless." We agree with this statement, as long as the 
components of T^ u are given in a coordinate basis, as 
e.g. T ^ in the proposal of Bicak et al [l3|. But we dis- 
agree with Einstein's objection, if the components are 
given in a LONB, T h y because these components can be 
directly measured by the FIDOs. The FIDO only needs 
clocks, meter sticks, and markers, which give the direc- 
tions of his spatial axes. The FIDO has only the metric 
of Special Relativity, r/mn = diag{— 1, +1, +1, +1}, avail- 
able in the tangent space at his space-time point. But 
the FIDO does not need any more structure, not a con- 
nection and not the metric potential functions g^^x), 
which in our case are given by the gravitomagnetic vec- 
tor potential A g , the solution and output of Einstein's 
equations. Therefore we disagree with Einstein's objec- 
tion quoted above (from which he concluded that one 
should no longer speak of Mach's Principle at all) , if the 
components of the energy-momentum tensor are given in 
a LONB. 



VII. THE LOCAL VORTICITY MEASURED BY 
NON-ROTATING OBSERVERS 

This quantity is defined as curlufl u ;d measured in the 
local inertial coordinate system which is comoving with 
the fluid, i.e. measured relative to the axes of local gy- 
rosopes ("local compass of inertia"). In general coordi- 
nates the local vorticity measured by non-rotating ob- 
servers is 

uo a = -e a ^ s u^ lUs , (38) 

where £ 0123 = — 1. Mach's principle, as formulated as a 
general hypothesis by Mach and made precise in Eq. 13211 , 
states that the gyroscope axes here follow the rotational 
flow J e of matter in the universe averaged with a r~ 2 
weight and an exponential cutoff at the H radius. In 
general the gyroscope axes here most definitely do not 
follow the motion of the local fluid here, i.e. relative 
to gyroscopes the local vorticity is nonzero according to 
Mach's principle in general. There is a special CcLSC, a 
rigid rotation of a fluid out to a perturbation radius 
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_R pGrt . In the limit Rpcrt/Rfi ~ * 00 the local vorticity 
measured by non-rotating observers vanishes. But this 
limit is uninteresting, since it produces an unperturbed 
FRW universe. — Unfortunately many authors have con- 
sidered the vanishing of the vorticity relative to the local 
compass of inertia to be a test for Mach's principle. See 



e.g. Ozsvath and Schucking's solution and discussion of a 
Bianchi IX model for perturbations of the original closed 
and static Einstein universe |Tfjj . In contrast we con- 
clude that the vanishing of the vorticity relative to the 
local compass of inertia is not relevant as a test of Mach's 
principle. 
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